Emergent multipolar spin correlations in a fluctuating spiral - 
The frustrated ferromagnetic 5 = 1/2 Heisenberg chain in a magnetic field 



O 

o 

O 

O 



i 

CO 



i 

£ 

O 

o 



(N 
> 

m 

On 

o 

oo 
O 

> 

X 



Julien Sudan, 1 Andreas Liischer, 1 and Andreas M. Lauchli 2 * 

1 Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), CH-1015 Lausanne, Switzerland 
Max Planck Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, D-01187 Dresden, Germany 

(Dated: October 9, 2009) 

We present the phase diagram of the frustrated ferromagnetic 5 = 1/2 Heisenberg J\ — Ji chain in a magnetic 
field, obtained by large scale exact diagonalizations and density matrix renormalization group simulations. A 
vector chirally ordered state, metamagnetic behavior and a sequence of spin-multipolar Luttinger liquid phases 
up to hexadecupolar kind are found. We provide numerical evidence for a locking mechanism, which can drive 
spiral states towards spin-multipolar phases, such as quadrupolar or octupolar phases. Our results also shed light 
on previously discovered spin-multipolar phases in two-dimensional S = 1/2 quantum magnets in a magnetic 
field. 

PACS numbers: 75.10Jm, 75.30.Kz, 75.40.Cx, 75.40.Mg 



Introduction Spiral or helical ground states are an old and 
well understood concept in classical magnetism [1], and sev- 
eral materials are successfully described based on theories of 
spiral states. For low spin and dimensionality however quan- 
tum fluctuations become important and might destabilize the 
spiral states. Given that spiral states generally arise due to 
competing interactions, fluctuations are expected to be partic- 
ularly strong. 

A prominent instability of spiral states is their intrinsic twist 
(Si x Sj) (vector chirality) [2], It has been recognized that fi- 
nite temperature [3] or quantum [4] fluctuations can disorder 
the spin moment (S,) of the spiral, while the twist remains 
finite. Such a state is called p— type spin nematic [5], In the 
context of quantum fluctuations such a scenario has been con- 
firmed recently in a ring-exchange model [6], while possible 
experimental evidence for the thermal scenario has been pre- 
sented in [7]. The twist also gained attention in multiferroics, 
since it couples directly to the ferroelectricity [8]. 

In this Letter we provide evidence for the existence of yet 
a different instability of spiral states towards spin-multipolar 
phases. The basic idea is that many spin-multipolar order pa- 
rameters are finite in the magnetically ordered spiral state, 
but that under a suitable amount of fluctuations the primary 
spin order is lost, while a spin-multipolar order parameter sur- 
vives. We demonstrate this mechanism based on the mag- 
netic field phase diagram of a prototypical model, the frus- 
trated 5=1/2 Heisenberg chain with ferromagnetic nearest- 
neighbor and antiferromagnetic next nearest-neighbor inter- 
actions. Furthermore we show that this instability provides 
a natural and unified understanding of previously discovered 
two-dimensional spin-multipolar phases [9, 10]. 

Hamiltonian and Methods To be specific, we determine 
numerically the phase diagram of the following Hamiltonian 




H 



•h 2_^ Si • Sj+i + Ji 2_^ Si • S,-| 



hJ2 S * 



(1) 



and we set Ji = — 1, J2 > in the following. Si are spin- 
1/2 operators at site i, while h denotes the uniform magnetic 
field. The magnetization is defined as m := 1/L £\ Sf ■ We 



FIG. 1 : (Color online) Phase diagram of the frustrated ferromagnetic 
chain (1) in the J1/J2 vs. m/m sat plane. The grey low-m region 
exhibits vector chiral long range order. The colored regions denote 
spin-multipolar Luttinger liquids of bound states of p = 2, 3, 4 spin 
flips. Close to saturation the dominant correlations are multipolar, 
while below the dashed crossover lines, the dominant correlations 
are of SDW(p)-type. The tiny cyan colored region corresponds to 
an incommensurate p = 2 phase. The white region denotes a meta- 
magnetic jump. Finally the scribbled region close to the transition 
J1/J2 — > —4 has not been studied here, but consists most likely of 
a low field vector chiral phase, followed by a metamagnetic region 
extending up to saturation magnetization. The inset shows the same 
diagram in the J1/J2 vs. h/h sa ,t plane. 



employ exact diagonalizations (ED) on systems sizes up to 
L = 64 sites complemented by density matrix renormaliza- 
tion group (DMRG) [11] simulations on open systems of max- 
imal length L = 384, retaining up to 800 basis states. 

Overall Phase Diagram The classical ground state of the 
Hamiltonian (1) is ferromagnetic for J2 < 1/4 and a spiral 
with pitch angle tp = arccos(l/4Ja) e [0, 7r/2] otherwise. 
The Lifshitz point is located at J 2 = 1/4. In a magnetic field 



the spins develop a uniform component along the field, while 
the pitch angle in the plane transverse to the field axis is unal- 
tered by the field. 

The zero field quantum mechanical phase diagram for S = 
1/2 is still unsettled. Field theoretical work predicts a fi- 
nite, but tiny gap accompanied by dimerization [12, 13] for 
J2 > 1/4, which present numerical approaches are unable to 
resolve. The classical Lifshitz point J2 = 1/4 is not renor- 
malized for S — 1/2, and the transition point manifests itself 
on finite systems as a level crossing between the ferromag- 
netic multiplet and an exactly known singlet state [14]. The 
theoretical phase diagram at finite field has recently received 
considerable attention [15, 16, 17], triggered by experiments 
on quasi one-dimensional cuprate helimagnets [18, 19, 20]. 
One of the most peculiar features of the finite size magneti- 
zation process is the appearance of elementary magnetization 
steps of AS Z — 2, 3, 4 in certain J2 and m regions. This has 
been attributed to the formation of bound states of spin flips, 
leading to dominant spin-multipolar correlations close to sat- 
uration. A detailed phase diagram is however still lacking. 

We present our numerical phase diagram in the J2/IJ1I 
vs. m/m sat plane in Fig. 1. At least five different phases 
are present. The low magnetization region consists of a sin- 
gle vector chiral phase (grey). Below the saturation magne- 
tization we confirm the presence of three different multipolar 
Luttinger liquid phases (red, green and blue). The red phase 
extends up to J2 — > 00 [16], and its lower border approaches 
m = + in that limit. All three multipolar liquids present a 
crossover as a function of m,/m sat , where the dominant corre- 
lations change from spin-multipolar close to saturation to spin 
density wave character at lower magnetization. One also ex- 
pects a tiny incommensurate p = 2 phase close to the p = 3 
phase [17], which we did not aim to localize in this study. Fi- 
nally the multipolar Luttinger liquids are separated from the 
vector chiral phase by a metamagnetic transition, which occu- 
pies a larger and larger fraction of m as J 2 — > 1/4+, leading 
to an absence of multipolar liquids composed of five or more 
spin flips. In the following we will characterize these phases 
in more detail and put forward an explanation for the occur- 
rence and locations of the spin-multipolar phases. 

Vector chiral phase For m > we reveal a contiguous 
phase sustaining long range vector chiral order [21] (breaking 
discrete parity symmetry), similar to phases recently discov- 
ered for Ji, J2 > [22, 23]. Direct evidence for the presence 
of this phase is obtained from measurements of the squared 
vector chiral order parameter: 



K 2 (r,d) := ([S x S d ] z [S r x S r+d Y 



(2) 



In Fig. 2 we display DMRG results for long distance corre- 
lations between J\ bonds (d = 1, black symbols) and J 2 
bonds (d = 2, red symbols) obtained on a L = 192 system. 
The three chosen values of J2 reflect positions underneath 
each of the three spin-multipolar Luttinger liquids shown in 
Fig. 1. The non-monotonic behavior of the correlations at 
very small m is probably a finite size artifact or convergence 
issue. Beyond the long range order in the vector chirality, 
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FIG. 2: (Color online) Squared vector chirality order parameter k 
(Eq. 2) in the low magnetization phase for different values of J2/I Ji| 
as a function of ra/m sat : (a) J2/I Ji|=0.27, (b) J2/I Ji|=0.32, (c) 
Jil\ Ji I = 0.38. The order parameter vanishes in the multipolar 
Luttinger liquids. 



the system behaves as a single channel Luttinger liquid (with 
central charge c = 1, confirmed by our DMRG based entan- 
glement entropy analysis [24]) with critical incommensurate 
transverse spin correlation functions [22]. The transition to 
the spin-multipolar phases at larger m seems to occur gener- 
ically via metamagnetic behavior, c.f. left and right panels of 
Fig. 2. For the parameter set in the middle panel we expect 
the same behavior, but it can't be resolved based on the sys- 
tem sizes used. 

Multipolar Luttinger liquid phases The Hamiltonian 
Eq. (1) presents unusual elementary step sizes AS Z > 1 in 
some extended J2/IJ1I and m domains, where AS Z is in- 
dependent of the system size [15]. This phenomenon has 
been explained based on the formation of bound states of 
p = AS Z magnons in the completely saturated state, and at 
finite 77i/m sat a description in terms of a single component 
Luttinger liquid of bound states has been put forward [16, 17]. 
We have determined the extension of the AS Z — 2, 3, 4 re- 
gions in Fig. 1, based on exact diagonalizations on systems 
sizes up to 32 sites and DMRG simulations on systems up to 
192 sites. The boundaries are in very good agreement with 
previous results [15] where available. The AS Z — 3 and 4 
domains form lobes which are widest at m sat and whose tips 
do not extend down to zero magnetization. The higher lobes 
are successively narrower in the J2 direction. We have also 
searched for AS" 2 = 5 and higher regions, but found them to 
be unstable against a direct metamagnetic transition from the 
vector chiral phase to full saturation. Individual bound states 
of p > 5 magnons do exist (see below), but they experience a 
too strong mutual attraction to be thermodynamically stable. 

An exciting property of the Luttinger liquids of p bound 
magnon states [17] is that the transverse spin correlations are 
exponentially decaying as a function of distance due to the 



binding, while p-multipolar spin correlations 
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11 0+n 11 S r +n) ~ ( 1) 



n=0 



n=0 



1/K 



(3) 



are critical with wave vector 7r (multipolar correlations with 
p' < p also decay exponentially), p = 2,3,4 corre- 
spond to quadrupolar, octupolar and hexadecupolar corre- 
lations respectively. Therefore they can be considered as 
one-dimensional analogues of spin multipolar ordered phases 
found in higher dimensions. Another important correlation 
function is the longitudinal spin correlator, which is also crit- 
ical [17]: 



(S S r ) 



(1 - m/m aa , t )Trr 



P 



K 



(4) 



We verified numerically c = 1 and determined the Luttinger 
parameter K as a function of to and J 2 by fits to the local 
S z profile in DMRG simulations. An important information 
is contained in the crossover line K = 1 where ^-multipolar 
and longitudinal spin correlations decay with the same expo- 
nent. This crossover line is shown for the three lobes in Fig. 1 . 
Close to saturation the spin-multipolar correlations dominate 
while towards the tip of the lobes the longitudinal spin cor- 
relations decay more slowly. The crossover line is rather fiat 
in the J2 direction, but drops to lower to values when going 
from p — 2 to 3 and 4. This can be understood from the 
analogy with a Luttinger liquid of hardcore bosons of bound 
states which effectively becomes more dilute when increasing 
p at constant m. 
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FIG. 3: (Color online) Wave vector q maximizing the longitudinal 
(inset: transverse) spin structure factor as a function of the magneti- 
zation m/m sa t for selected J2/I Ji| parameters, each one belonging 
to a different p sector at high m. The straight lines are guides to the 
eye and show n(l — m/m sat )/p. The different symbols denotes ED 
results on system sizes ranging from 18 to 28 sites. 

Interpretation After having established and characterized 
the phase diagram, we are lead to the intriguing question 
whether there is a deeper relation between the spin multipolar 
phases which would enable us to understand their presence, 
and possibly predict related phenomena in other systems. 



Let us first investigate how the longitudinal and transverse 
equal-time spin structure factors S zz (q) and S xx (q) evolve 
as a function of to. In Fig. 3 we display the location of the 
maximum of S zz (q) and S xx (q) (disregarding the q = 
peak in S zz {q) due to the total magnetization) for three rep- 
resentative J2 values. At m = it is known that S(q) has a 
maximum at an incommensurate position g max (</2), which is 
strongly renormalized compared to the classical expectation 
9max S (^2) = arccos(l/4J 2 ) [25]. In the low magnetization 
region, corresponding to the vector chiral phase, the location 
of the maxima of both structure factors are only weakly de- 
pendent on m, and in a first approximation remain the same 
as for to = 0. However as to is increased, the g max of S zz (q) 
locks onto a straight line with slope — ir/p. It seems that if 
9max(^2)/ 7r > 1/3 at m = then the magnetization process 
will enter the p = 2 multipolar phase at larger m. If instead 
1/3 > g , max(^2)/7r > 1/4 at to = then the system will en- 
ter the p — 3 phase. Based on an extended analysis including 
several more J2/J1 values we are thus lead to conjecture that 
if 



1/P > 9ma X (J 2 )/7I- > I/O +1) at TO = 0, 



(5) 



then the higher to region will lock onto the line with slope 
— it /p. According to the behavior of the longitudinal spin cor- 
relations in the multipolar Luttinger liquids (Eq. 4), the multi- 
polar liquid of order p leads precisely to a slope of —n/p. 

A complementary point of view is provided by the consider- 
ation of a simple classical spiral. The corresponding magnetic 
structure factor has a peak at the propagation vector q of the 
spiral. It is straightforward to show that in such a state the in- 
plane multipolar correlations of order p (analog to Eq. 3) have 
a propagation vector q p = p-q. Of course in a magnetically or- 
dered state these multipolar correlations are not primary order 
parameters, but strong quantum fluctuations can wipe out the 
spin order, while a particular multipolar order parameter sur- 
vives. Applying this observation to the structure factor results, 
we can understand the rule (5) as a locking of the multipolar- 
p correlations to ir, which precisely corresponds to the stag- 
gered nature of the multipolar-p correlations in Eq. 3. 

Despite the fact that there are no stable multipolar Luttinger 
liquids with p > 4, it is possible to examine the process of 
individual bound state formation by pushing previous calcu- 
lations from p < 4 in Ref. [17] to p = 8. We performed 
ED without truncation on system sizes up to L = 64, cal- 
culate the binding energy of p flipped spins [expressed as 
h sa t(p) I h sa t(p = 1)] and display it in the upper panel of 
Fig. 4. We obtain a surprising series of stable bound states 
up to p = 8, which are successively smaller in their J2 ex- 
tent, suggesting that J 2 = 1/4 is an accumulation point of 
p — > 00 bound states. We now proceed to a comparison of 
the stability domain of the bound states to those predicted by 
the locking rule (5) applied to the zero field structure factor. 
We fit a power law g max (J 2 ) ~ ( J 2 - 1/4) 7 (with 7 sw 0.29) 
to the finite size transition points, as shown by the bold black 
line in the lower panel of Fig. 4. This fitted line is used to con- 
struct the boundaries, giving rise to the lower color bar. The 
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FIG. 4: (Color online) Upper panel: ED results for the domain of 
stability of bound states of p spin flips above the saturated ferro- 
magnetic state. Lower panel: Location g ma x of the maximum of the 
zero field spin structure factor 5(g) as a function of J2. The colored 
lines indicate the construction of the domain boundaries, based on 
the locking rule 
classical model. 



agreement between the transition boundaries in the upper and 
lower panels is excellent (apart from the p = 2 to p = 3 tran- 
sition). It is striking that we are able to reproduce the domain 
of stability of bound states of magnons at the saturation mag- 
netization, based solely on the spin structure factor obtained at 
zero field. This constitutes strong evidence for an approximate 
validity of the locking rule Eq. 5 and the presence of a lock- 
ing mechanism which pins the multipolar correlations to n at 
larger m as opposed to incommensurate transverse spin cor- 
relations for lower m. The transition from p — 2 to p = 3 is 
somewhat shifted compared to the prediction, and we attribute 
this discrepancy to the formation of a peculiar incommensu- 
rate p = 2 bound state [17, 26], leading to a collapse of the 
kinetic energy of the p — 2 bound state. 

Our considerations presented above are not restricted to one 
dimension. Indeed the recently reported quadrupolar [9] and 
octupolar phases [10] in frustrated 5=1/2 ferromagnets 
in a field also allow an interpretation as fluctuation destabi- 
lized magnetic states. In the J\-J<z model on the square lattice 
there is a magnetically ordered striped phase [Q = (71", 0)] 
in the neighborhood of the "bond nematic" phase. The near- 
est neighbor bond quadrupolar correlation functions in such 
a classical striped state have momentum (0, 0) and are odd 
under a 7r/2 rotation of one bond. Interestingly this is ex- 
actly the signature of the correlation function reported in the 
magnetically disordered, "bond nematic" phase. Similarly for 
the "triatic" or octupolar phase, there is a canted magnetic 
three-sublattice state [Q = (47r/3, 0)] neighboring the "tri- 
atic" phase at high fields [10]. The derivation of the octupo- 
lar correlations on a triangle in the classical state again yields 
octupolar correlations matching the symmetry of the bound 
states reported in the "triatic" phase. Based on the success of 
our considerations we speculate that it could be possible to 



stabilize a uniform hexadecupolar phase in the square lattice 
Ji — J3 model with ferromagnetic nearest neighbor interac- 
tions. 

Experimental relevance The different phases we found 
should be detectable experimentally in the quasi- ID com- 
pounds studied in Refs. [18, 19, 20]. In neutron scattering 
experiments these phases are expected to manifest as follows: 
i) a Bragg peak at low magnetizations corresponding to spi- 
ral order in the plane transverse to the magnetic field, ii) a 
Bragg peak at intermediate magnetizations at wave vector 
7r(l — rn/m sa t)/p corresponding to SDW(p) magnetic order 
modulated along the field and iii) the absence of nontrivial 
magnetic Bragg peaks at high fields in the spin multipolar or- 
dered phases. 

Conclusions We have established the phase diagram of 
the frustrated ferromagnetic S — 1/2 Heisenberg chain in a 
uniform magnetic field. We provided an explanation for the 
appearance of a large number of bound states upon approach- 
ing the Lifshitz point at J2 = 1/4, based on a locking mech- 
anism leading to spin multipolar phases. It is an open prob- 
lem whether these fluctuation driven multipolar phases and 
the locking mechanism also appear for S > 1/2. In a recent 
paper by Hikihara et al. [27] a phase diagram very similar to 
ours has been obtained. 
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